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Abstract. The many-body theory of photoemission in solids is reviewed with emphasis on 
methods based on response theory. The classification of diagrams into loss and no-loss diagrams 
is discussed and related to Keldysh path-ordering book-keeping. Some new results on energy 
losses in valence-electron photoemission from free-electron-like metal surfaces are presented. A 
way to group diagrams is presented in which spectral intensities acquire a Golden-Rule-like form 
which guarantees positiveness. This way of regrouping should be useful also in other problems 
involving spectral intensities, such as the problem of improving the one-electron spectral function 
away from the quasiparticle peak. 



1. Introduction 

Photoemission spectroscopy ^ (PES) is since long a major tool for gaining information of 
the electronic structure of matter. For instance, angular resolved photoemission spectroscopy 
(ARPES) makes it possible to measure the occupied quasiparticle bands, and photoemission from 
core levels yields important local information at the atom being excited. Owing to the usually 
rather short mean free path, PES is also well suited for studying surface states, adsorbates at 
surfaces, etc. 

Photoelectron spectra are usually interpreted in terms of the one-electron spectral function 
A(e) corresponding to a sudden removal of an electron. In emission from solids, however, the 
photoelectron is created at a certain distance from the surface and may undergo losses on its way 
out. In the sixties, Berglund and Spicer|2j proposed a semi-empirical model involving primary 
excitation, transport to the surface, and transmission through the surfaces as three separate 
steps. Clearly a three-step description can at best be an approximation. An important step 
forward was taken by Schaich and AshcroftjH] and by Mahan[3], who formulated the problem as 
a one-step quantum-mechanical process and presented model results for independent electrons. 
The Schaich-Ashcroft response formulation was soon generalized to account for interaction by 
Caroli et aZ.jH] who applied the formalism to phonon effects. Among important early works I 
would also like to mention the studies of plasmon losses by Langreth and collaborators [HI 0012] • 
In particular, they explain why the satellites are greatly modified and suppressed also when the 
excitation energy is quite high, and why the satellites are so weak in, say, x-ray emission. 

In strict one-electron theory, the photoelectrons cannot undergo losses by scattering against 
other electrons, phonons, and impurities. As a consequence, the entire solid will contribute to the 
photocurrent which is then limited by the penetration depth of the electromagnetic radiation. 
In reality, the current is limited by the escape depth of the photoelectrons, in practice usually 
just a few atomic layers. Thus the mean free path must somehow enter the theory, and in the 
mid-seventies ways to remedy this shortcoming were developed by Feibelman and Eastman 10 , 
Pendrv[TT]. and Liebsch fT!?! I13| . In essence, the no-loss photocurrent is to be calculated in 



a one-step but independent-electron way with one important exception: The photoelectron is 
assumed to move in the non- Hermit i an optical potential, the (advanced) self-energy S . The 
photoelectron orbitals then get damped inside the sample, and the yield becomes limited by the 
photoelectron escape depth as it should. Today there exists several computer codes based on 
these ideas which take the underlying bandstructure into account in a realistic way. 

The crucial point of introducing damped photoelectron orbitals was first obtained by 
LangrethjHj in the case of impurity scattering, and by Caroli et a/. [21 in the case when the 
escape depth is limited by phonons. In the mid-eighties, by Almbladh|14j and Bardyszewski and 
Hedin ^S] considered the problem of interacting electrons using two seemingly rather different 
techniques, one based on many-body perturbation theory and one based on scattering theory. 
These works largely justify the above one-electron-like methods for the no-loss photocurrent but 
they also predict corrections, which, for instance, modify the optical matrix elements. 

There are several excellent recent reviews of solid-state photoemission |16l 1171 ITS] . In this 
paper, I would like to concentrate on how the techniques of non-equilibrium Green's functions 
can be used both for obtaining approximations suitable for calculations in real materials as 
well as for obtaining general results like the structure of the quasiparticle no-loss part of a 
photoelectron spectrum. 

2. Basic theory 

I begin by a brief account of the quantum-mechanical description of the steady photoelectron 
current. As mentioned above, methods based on scattering theory ^3HJ and methods based 
on quadratic response formalism and many-body perturbation theory [3} 15] I20| 16] 19] I14j . I begin 
by one-particle theory in order to illustrate the various limiting procedures, and generalize then 
to interacting systems. 

2.1. Independent particles 

Let us first consider independent electrons. The photoemission process can be considered as 
an inelastic scattering process with a photon in the initial state and an asymptotically free 
electron in the final state. In order to strengthen the analogy with scattering, let us take the 
system to occupy an finite volume v enclosed in a quantization box of volume V. Because the 
photoelectron enters in the final rather than the initial state, its wavefunctions will obey ingoing 
wave boundary conditions |2 1 \ HH] . To first order in the radiative coupling we obtain 



Here, J p is the momentum-decomposed photocurrent, A the coupling to the radiation field (cf. 
Eq. LP the photon energy, and e u one-electron energies. Atomic (Hartree) units are used. 

The photoelectron spectrum probes the occupied part of the one-electron spectrum. For 
independent electrons, the spectral function can be written 
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in terms of eigen-orbitals cp n . In terms of the spectral function, we have 
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for photoelectron energies below threshold n+u. Here \x is the Fermi level relative to the vacuum 
level. 



In order to compare with the response formulation of Schaich and Ashcroft, let us introduce 
orbitals 4> p which are plane waves well outside the sample and which vanish inside, and let c p be 
the associated destruction operator. We assume that the external field has been slowly switched 
on in the remote past, (modulated by e vt ), and take the average 2r)(cpC p ) as a measure of the 
steady momentum-decomposed photocurrent. The details how the plane waves are truncated 
inside and near the sample will not influence the final results. We have 

K(t) = H+ \5He~ iujt + 5H^e iujt ] e vt , (4) 

where H describes the system in the absence of external fields, and where 

6H = - /"j(r)-A(r)=5>+ Q A w (5) 
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gives the coupling to the radiation field. We use a radiation gauge such that the scalar potential 
5& vanishes. (The diamagnetic coupling then gives no contribution to leading order, see Ref [3].) 
We assume that the system has been subject to the field for t < 0, and define the photocurrent 
by 

J p = lim . 2r]{c±cp). (6) 

As in scattering theory the infinite volume limit should be taken before the limit ij — > , i.e., 
77 should be considered large compared to the level spacing due to the quantisation box V. 
Following Liebsch|T2j, we develop each initial-state orbital n in time to obtain 

<t>n, v ( t = °)= l + G r (e n + uj + ir])A + G r (e n -u + ir))A^ <j> n , (7) 

where G r ^ a (e) = l/(e — h±i5) are the retarded/advanced one-electron Green's functions, and h^i 
the one-electron Hamiltonian matrix of the unperturbed system. We next form the expectation 
value (cpCp) and sum over occupied states initial states to obtain 
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Here, G< is the usual lesser function, <A;|G r< (t) |Z) = i(N\c^ (0)ci(t)\N) , and \N) the iV-electron 
ground state. In order to obtain Eqs. ((HI El) we have used the fact that electrons in initial 
states n are confined to the solid, and we have dropped the negative-frequency term which is 
easily seen not to contribute to the photocurrent in the limit 77 — > + . In this limit one can 
further replace the truncated orbital (ftp with a plane wave <j) p (see e.g Ref. 14 J. Eq. Q is then 
equivalent to the Schaich- Ashcroft independent-electron result. 

In order to see the equivalence between the scattering and response formulations, we express 
Qr/a - m t erms Q f Green's function G r J a for free space and the scattering matrix T r / a , 

G r ' a {e) = G r /a (e) + G r J a {e)T r > a {e)G T J a {e) . (10) 

As remarked above, we can replace the truncated plane wave <j)p with <p p in the limit V — > 00, 
77 — > + . In this limit we have 

(n|l + AG a {e n + oj -irf)\p) = (n\AG a (e n + lj - irj)\p) 

= — — -(n\A[l + G a (e p -iri)T a (e p -i V )]\p) 

e n +uj - e p -vq 



and the wave operator 1 + GgT a transforms 4> p into x p ■ I n this wa y> 
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| (p\G r (e n + w + zr/)A|n)| 2 -» ' Up J " 2 - 2tt| ( Xp | A|n) \H((e n +u-e p ), 

which establishes the equivalence between the response and the scattering approaches for 
independent electrons. 

The geometry with a finite but possibly macroscopic sample is convenient for performing 
the somewhat subtle limits above, but in actual calculations one usually let the sample grow 
until it fills a half-space. The final-state orbitals x P are then replaced by their two-dimensional 
(2D) analogues, i.e., time-reversed states for scattering against the surface of the sample. These 
states are usually referred to as time-reversed LEED states because similar 2D scattering enters 
in the problem of low energy electron diffraction (LEED). 

As stressed in the Introduction, strict one-particle theory is unphysical and leads to a yield 
limited by the radiation penetration depth rather than the photoelectron escape depth. If the 
propagators in Eq. © are replaced by their interacting counterparts, the final-state orbitals get 
damped inside the sample and the photocurrent limited by the escape depth as it should. 



2.2. Interacting electrons 

We expand (cpC p ) directly in orders of the perturbation. The linear response vanishes identically 
because all electrons are confined to the sample in the ground state, c p c p \N) = 0. In the 
radiation gauge, the diamagnetic coupling involving A 2 gives no contribution. The remaining 
contributions can be written 

pO p0 

J p {e) = lim 2r, / dt dt'e^*^ e^ t+t>) (N\5H^ (t)N p (0)5H(t')\N) . (11) 
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Eq. (jlll) involves a three-particle path-ordered Green's function 

R P (t,t') = Y, A^ n A Mn (iV| c +(0 C „(t)ct(0)c j5 (0) C +(t')c m (t , )|A r }, (12) 
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which may be evaluated by standard Keldysh technique |22j. To lowest order, all unconnected 
contractions vanish because 

G<(t)\p) = 0, (p\G<(t) = (13) 

both with and without interaction. By the same argument, there is only one non- vanishing 
connected contraction (see Fig. Q), which gives 

4°>(f,0 = i(p\G>(-t')AG<(t-t')A^G>(t)\p) 

= i(p\G r (-t')AG<(t-t')^G a (t)\P) (14) 

when t, t' < 0. Eq. (|14|) gives back the previous independent-electron result in Eq. Q. 

In the present problem, the correlation R p (t, t') in Eq. ()12|) is a ground-state correlation which 
simplifies the book-keeping of Keldysh diagrams. The Keldysh propagators refer to the contour 
in Fig. @ which means that we assume that the interacting ground state can be generated by 
switching on the interaction adiabatically in the remote past, 



R p (t,t') = (N,0\S^ (oo)5H^ (t)ct S(oo)S^ (oo)c p 5H(t')S(oo)\N,0) . 
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Figure 1. The lowest-order non- vanishing diagram for the photoemission process 



Thus t = t-, t' = t' + , and CpCp is straddling the two branches of the Keldysh contour in Fig. 
@. (The indices ± refer to the positive and negative time-ordered parts, respectively.) The 
propagators joining points fi, v of the different time-ordered parts are 

G c (t-t'), fi, v = + 
, , _ G c (t-t'), fx, u = - 

^[t^tu)- Q<(j._ t >) } M=+)I/ = _ 

G>(t-t'), n=-, v = + , 

where G c , G c are propagators with positive and negative time-ordering, respectively, and G> 
the usual lesser /greater functions. At zero temperature we have 

a^(e)6^(e) = (16) 

for any Fermi propagators a and b. Further, with the above assumption of adiabatic switching, 
the expansion of a time-ordered propagator will only involve propagators with the same time 
ordering. 
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Figure 2. Keldysh contour at T = 



The contribution from a general diagram can be reduced to a form analogous to the one- 
electron Golden Rule form in Eq. In the expansion in full Green's functions, any diagram 
ends with a pair of electron lines joined to the vertex N„ representing the analyser (Fig. 
This means that the correlation R p in Eq. (|12j) is of the general form 

R p (t,t') = -J dxdrj dx' dT'{p\G c (-T')\x')(x'\I(t,t';T,T')\x)(x\G £ (T)\p), (17) 

where /, the remain when the exit fermion lines have been removed, is invariant under a 
simultaneous translation of all four times and can be written 

lit, t'; r, t') = f dvdV ' d f v'- v\ v> + v" - v) e -ivt+nSt> e i 1 />T-(vW-v")T\ (18) 
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We notice that G c and — G c can be replaced by retarded and advanced propagators, respectively, 
since p projects out the particle part. The action of G r l a on the orbitals p and the limit r\ — > 
can now be done analytically in much the same way as for independent particle, which gives 

J P = (Xp\I{u,u;e p: ep)\Xp)- (19) 

Here Xp is an one-electron time-reversed LEED orbital solved in the potential Vc + E a (e p ). (Vc 
is the total Coulomb potential from the nuclei and the equilibrium electron density.) The self- 
energy S a (e) is non-Hermitian with a finite imaginary part inside the sample. This makes Xp 
decaying inside the sample with a decay length determined by energy losses. Thus all parts of 
the spectrum including the no-loss part is now limited by the escape depth of the photoelectrons. 




Figure 3. Reduction of a general photoemission diagram 



In the case of optical responses and transport, the ideas of self-consistent approximations 
put forward by Kadanoff and Baym 23J [2IJ US] has been extremely successful. It is the 
present author's opinion that this level in not possible to achieve in the case of photoemission. 
In particular, the Bethe-Salpeter equation which emerges as a result of selfconsistency will 
include interaction lines of all orders. If we open up such a diagram corresponding to a given 
approximation to the self-energy, infinitely many photoemission diagrams will result. For sp- 
bonded solids, the most important conserving approximation is the self-consistent screened 
exchange approximation [251 I26| . usually referred to as the GW approximation. Fortunately, 
the fulfilment of macroscopic conservation laws are less crucial for useful approximations to 
photoelectron spectra. In photoemission, only a small fraction of the emitted photoelectrons 
are actually measured. The majority of photoelectrons which make up the overall background 
is usually of no interest at all; instead the interest is focused on the spectral shapes in a limited 
energy window. Often, processes of different orders will contribute at different energies. Losses 
to plasmons is here a typical example. What is important is to have approximations which can 
explain the observed spectral profiles and which fulfil basic requirements such as positiveness. 

2.3. Low-order diagrams 

Let us expand the basic 3-particle correlation R p in orders of the screened interaction W{uo) = 
fe _1 (Lj). The non- vanishing first-order diagrams are shown in Fig. ©. As was the case 
for independent electrons, the condition in Eq. (|13|) limits the number of non-vanishing 
contractions. Vertices belonging to the forward (backward) time-ordered parts are are labelled + 
(-), respectively. The diagrams (a-c) involve only self-energy insertions and should be omitted if 
we expand in the full Green's function. The remaining first-order diagrams fall into two classes, 
namely those with interaction line joining point in the same leg (d and e) and those with the 
interaction line joining points in different legs (f, g, and h). Consider first diagram f. It is well 
known that this diagram describes an extrinsic loss, but to illustrate the technique we give some 
details of the analysis. The interaction line here represents the boson-like correlation 

/•oo 

W>(rt,rV) = -i(N\5VH(rt)5V H (r't')\N) = -i / du/B(r, r', w / )e~ <w ' ( *~*' ) > (20) 
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where the spectral function B is given by the dynamical structure factor S(t,t',uj) convoluted 
with two Coulomb interactions (v), B(u) = v*S(lo) *v. The fermion line from t' to t represents 
G < . and the remaining lines time-ordered or anti-time-ordered functions G c l c . After reduction 
according to Eq. (fT9]> we obtain 

/0 roc 
de / dJ8(e + oj-J -e p ){x p \G c {e + uj)AA{e)A^G c \e + uj)\ Xp ) (21) 
-oo JO s V ' 
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where we have used an underbrace to indicate how the spatial coordinates in B{uj) should 
be connected to the electron propagators. Here and in the following we use the identity 
G c (e) = —G c ^(e) in order to express final results in spectral or lesser functions and functions 
with positive time-ordering. 

The obvious physical interpretation of Eq. (|21|) is that it describes primary photo-excitation 
from an initial occupied level (e) to a level of energy e + u. The photoelectron then suffers an 
energy loss by exciting a plasmon or electron-hole pair of energy u/ and hits the detector with 
an energy e + lo — u'. Consider next the diagrams g and h in Fig. 3. Using a similar technique 
as above we obtain 

de / dJ5(e+uj-J -e p )( Xp \G c {e + u)AA{e)G c \e-J)ti\x p ) (22) 

-oo JO s V ' 

describing loss processes where interference between the potential from the hole left behind and 
the outgoing photoelectron interfere El E] The 'intrinsic' losses, finally, are describe by the 
satellite structure of the hole spectral function A. The diagrams f-h plus the first-order satellite in 
A (c) give the first-order contribution to the first plasmon satellite as obtained by Inglesfield 27 
in the case of photoemission from core levels. 

We now turn to the first-order diagrams (d) and (e) with no interaction lines joining the (+) 
and (-) parts. In this case the interaction lines represent time-ordered correlations, 

/0 ("oo 
de / dLo'5(e + u;-e p )(x p \A^(e + u;,e)A(e)^\Xp) (23) 
-oo Jo 

where 

/T / 
— (n| G c (e - J)AG c {e' - u') \m) (24) 
2tt „ ' 

W c {w') 

is a first-order correction to the usual (time-ordered) vertex function. The diagram (e) gives 

the complex conjugate [J p d ^]* . The appearance of 5{e + uj — e p ) in Eq. ((2*3*)) shows that they 
contributes to the no- loss current. 

The remaining two diagrams (k) and (/), finally, describe a screening of the photon field in 
the solid [TO]. The screening will modify the electromagnetic field by local-field effects which are 
lattice-periodic deep inside the solid but more complicated near the surface. Most photoemission 
calculations to date do not take these effects into account. Effects of such screening has been 
considered by Liebsch 28 , Schattke,2""Q [SU], and by others. In the following discussion they will 
be considered included in the optical vertex A. 

With independent-electron propagators, the photoelectron orbital is approximated to zeroth 
order in the interaction and is undamped inside the solid. If we regard the diagrams skeletons 
and use full Green's functions (diagram (c) should then be omitted), the photoelectron orbital 
gets damped inside as it should. 
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Figure 4. First-order diagrams for the photoemission process 



2-4- Treatment based on scattering theory 

Scattering theory furnishes an equivalent description of the photoemission process and has been 
the basis for several interesting new results [HD E21- The obvious generalisation of the one- 
electron result is to replace the one-electron scattering orbital \xp~) by its many-body counterpart 
\f) = \Xp ; s ~) = l-^V — l,s, p,—) of energy Ef = e p + E S (N — 1) with an asymptotically free 
electron p and a N — 1-electron system left behind in some excited state s, 



j p = 2nj2\(f\sH\m\ 2 m-E f )) 



(25) 



The equivalence with the response formulation can be demonstrated with a variety of methods. 
One way is to consider the limiting behaviour of the amplitude 
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inherent in the response formulation and verify that it tends to 

(f\5H\N) 



Ei-E f + in 

when n — > 0. There are some subtleties involved because the photoelectron is identical the 
electrons in the sample. For a discussion, see e.g. Ref. |14j . 

A drawback with the above treatment is that every possible final many-body state has to be 
treated explicitly, a clearly impossible task in a 'full' treatment. In practice, however, one has 
so far only been able to account for losses to low orders. In sp-bonded materials, for example, 
the important excitations are boson-like and consist of plasmons and particle-hole pairs. The 
plasmons carry the major oscillator strength. Going back to Fig. (|4|). diagrams (c) and (/ — h) 



all have a real boson in the final state. By summing over possible intermediate virtual states 
with some subspace, the propagators become renormalised and one obtains an approximation 
to the first boson loss satellite with damping properly included. The idea to work in subspaces 
is the essence of the Hedin-Bardyszewski approach |15j. To make the problem more clear-cut, 
Hedin following an idea of McMullen et al. 33 considered the photoelectron as distinguishable 
from other electrons in the sample, the 'blue-electron model'. 

The very existence of the Golden-Rule expression in Eq. (|25|) suggests that also the results 
from response theory can be put in a similar form. This is indeed the case, and a technique how 
it may be achieved will be discussed at the end of this paper. 

3. No-loss current 

Photoemission spectra typically consist of a quasi particle-like peak and various losses at lower 
kinetic energy. A mathematically precise definition of a no-loss part can only be given in the 
uninteresting case where the sample ends up in the N — 1-particle ground state. If we are willing 
to accept a certain fuzziness, we can define a "no- loss" part with about the same precision as 
quasi particles below the Fermi level. In practice there is usually no difficulty to identify, say, 
core-electron quasi particles, or valence-electron quasi particles well below the Fermi energy so 
long as their lifetime broadening is small on the energy scale of interest. 

Our foregoing analysis suggest that interaction lines which join the positive and negative 
time-ordered parts of a diagram correspond to energy losses in the final state. There must 
always be at least one G < line joining the two optical vertices in a diagram, i.e, there is always 
an excitation where an electron has been removed. In first order, the diagrams (c, / — h) in Fig. 
(jlj) contribute to additional losses involving particle-hole pairs or plasmons because they involve 
W > . In an interaction-picture representation for 3-current correlation R p we have 

R p (t, t') = (0^ [S(-oo, oo)5H(t)c{0)] T [S(-oo, oo)c(0)SH(t)] |0>. (26) 

We insert a complete set of final N — 1-electron states between the forward- and backward time- 
ordered parts and expand. The N — 1 electron ground state corresponds to ejecting an electron 
from the Fermi surface adiabatically. If the sample has well-defined quasi particles it is natural 
to filter out all final states with energies a couple of lifetime widths away for the quasi-particle 
energy. This corresponds essentially to excluding interaction and electron lines which cross the 
two time-ordering parts. If we sum all diagrams of this kind, we obtain 

J p = 2vr ^ 5(e p -e-u;)(x-\Me + L0,<i)A(e)A\e + u;,e)\x-) de (27) 

J — oo 

(see Fig ©). Here, A(e,e') is the screened, time-ordered vector-coupling vertex. In order to 
obtain only the no-loss part the integration with respect to e should be confined to the quasi 
particle part of A(e). How this should be done cannot be specified exactly, as a consequence 
of the fact that the sample is generally left in an unstable excited state. The impreciseness 
is of the order of the lifetime width of the hole left behind. Experimental angular-resolved 
spectra show that the quasi-particle concept is useful also for states far below the Fermi energy 
and that spectra of real systems can be considered composed of a no-loss part that reflects the 
quasi-particle bandstructure, and a loss part that reflects the dynamical response of the system. 

The vertex correction is clearly related to the problem of how the optical matrix elements 
should be calculated. In the dipole approximation, the velocity (A&; ~ {k\A. ■ p\l}), length 
(Afc; ~ (fc|E • r\l)), and acceleration (A^ ~ (fc|^4 • Vw(r)\l)) formulas are equivalent. For 
independent particles the equivalence readily follows by commuting the one-electron Hamiltonian 
with the momentum, and w is then the effective one-electron lattice potential. When the 
independent-electron propagators are replaced by their interacting counterparts, however, the 



equivalence is lost. There is also a question how the effective potential w should be chosen. 
This problem resolved by the present author in the mid-eighties |34 j . The replacement of 
the bare optical matrix element by the dressed vertex A v /x,/ a restores the equivalence almost 
completelv|34). The remaining discrepancy is of order of the lifetime width of the hole left 
behind and is clearly related to the fuzziness of quasi particles away from the Fermi surface. 
This again shows that Eq. (|27|) furnishes a good description in systems with well-defined quasi 
particles. The acceleration matrix elements involve e _1 (w)A • Vw nuc (r) (a polarisation where 
are nuclei have been shifted a distance A. The last result was first obtained by Hermeking 20 . 

A 



N p 



Figure 5. The no-loss part of the photocurrent 



4. Energy Losses 

We shall here mainly be concerned with plasmon losses in sp-bonded materials. As is well known 
plasmon losses can be produced either "intrinsically" in connection with the primary photo- 
excitation, or " extrinsically" by the escaping photoelectron on its way out. At lower kinetic 
energies the two modes of plasmon excitation interfere, leading to a suppression of the satellite 
strength. This interference has in the past been described using semi-classical models where the 
photoelectron is assumed to follow a classical trajectory and where the photoelectron and the 
hole left behind are represented by external potentials EHl EZ1 EE]- I n the case of core-electron 
photoemission, Hedin et al. [611 have shown that a fully quantum-mechanical description of 
the leading satellites agrees quite well with a semi-classical description for excitation energies 
a couple of plasmon energies above threshold. It also shows that interference between intrinsic 
and extrinsic plasmon production persists up to the keV range, thereby confirming the earlier 
works by Chang and Langreth[3 |Hj. However, there are to our knowledge no similar studies for 
photoemission from valence states. 

As will be discussed below, a straight-forward expansion in orders of the screened interaction 
W will not in general guarantee positive spectral densities. I will here outline a different 
summation order in which the partial sums have a golden-rule-like form involving transition 
amplitudes r/j. Thus, each partial sum has the generic form 

I(oj) = Y,\Tfi\ 2 5(E i +u J -I f ), (28) 
/ 

where i, f label real excitations in the initial and final states, respectively. The transition 
amplitudes can be related to cut Keldysh diagram. By performing an expansion of the 
amplitudes to a given order, each partial sum gives a positive result. The procedure amounts 
to splitting complete Keldysh diagrams in pieces and to sum them a particular order such that 
every partial sum is positive. If carried out of infinite order, 'all' complete diagrams are summed. 
I will outline the procedure with the low-order plasmon losses in photoemission as an example, 
and give some new model results for valence photoemission from a simple-metal surface. 




4-.1. Lowest- order contributions 

Let us begin with the lowest-order contributions where no real losses have occurred. Up to first 
order in W, the zeroth-order diagram in Fig. Q an d diagrams (d) and (e) in Fig. 0] contribute 
and will involve the combination 



in which the transition amplitudes rather than the spectral intensity has been truncated to a 
given order. 

Contributions with no real excitations in the final state are somewhat artificial and will give 
no contribution when particle-hole excitations are included in the screened interaction. The 
mechanism is essentially the same as in the Langreth theorem j4Uj. In the case of the simpler 
G < function we have 



If there are no real excitations in a particular spectral region, S < vanishes, and the first term in 
Eq. (|29[) gives the entire contribution. As soon as real excitations occur, the first term tends to 
zero in the infinite-volume limit and the entire contribution comes from the second term. What 
actually occurs is that previously sharp single-particle- like excitations become broadened. By the 
same mechanism, contributions with no real excitations in the final state except a perfectly sharp 
quasi particle will not contribute to the photoemission current when broadening by particle-hole 
pairs is included. 

We next proceed to the first-order diagrams (c) and (/ — h) in Fig @ which contribute to 
a one-boson loss. These four diagrams nicely combine into a complete square of intrinsic and 
extrinsic satellite amplitudes. We make use of the fact that the spectral functions B(uj) and 
A(e) are Hermitian and positive-definite. They can thus be diagonalised, 




( X ;\[A + AW}A(e)[A + A^\ Xp ) 



G<(e) = [l + G r (e)E r (e)]G, 



(e)[l + £ a (e)G a (e)] + G r (e)E<(e)G a (e). 



(29) 




(30) 




(31) 



The "intrinsic' part (c) can be written 



4 c )=27r< X -|AA<(e p -u,)At| x -), 



(32) 



where A < = G C T < G C ^ and T < = — iT, < /2-k. In the GW approximation j2o] ■ we have 
(r|r<( ep -a,)|r'} = 



^ / duj'al(e p +J- w)^(u/)u n (e p + u)' -u, r)u* n (e p + J - r')u M (a/, r)v*(uj', r) 



We insert this in Eq. (|21|) to obtain 




(33) 



where 



M»£ a = (x;|AG c (e p - u;)\v;(u')u n (e p + J - u)) . 



(34) 



We next turn to the purely extrinsic part (/). Again we rewrite A < and B{u) in a diagonalising 
basis to obtain 

poo 

J^=2nJ2 ^|M^0lMfeW- W )&>'), (35) 



n/A 



where in this case 



1,1,6 



(x p v fl (u')\G c (e p + u')A\u n (e p + J - u)). 



(36) 



The two interference terms (g — h), finally, can be expressed in a similar way but with matrix 
elements a M\ j if ) . The total contribution can then be written 



n/A 



duo' 



ilffiu/) a 2 n {e p + J - u;)b 2 m (u') 



(37) 



where 

M M = Mi,i, a + M lxb . (38) 

The approximation of keeping only the renormalised diagrams (c, / — h) has a close 
correspondence to the GW approximation and may be termed a l GW approximation for 
photoemission'. 

In the simplest non-selfconsistent GqW approximation, the spectral function in G < should 
be approximated by its independent-particle equivalent in S^. (For different levels of 
selfconsistency, see Ref. 030 I n this approximation, a^ n {e) is replaced by 5(e — e m ) in the above 
expressions, and remaining propagators by their GqW equivalents. One can also in principle 
do partial selfconsistency with respect to G, which in this context amounts to calculate A(e) 
selfconsistently from Eq. (|31jl. the resulting £>, and the Langreth theorem, Eq. ([29 j). 

The matrix elements M\ t \ a and Mm can be represented graphically at cut diagrams with 
only + vertices and one dangling interaction line W > as in Fig. |f3 When joined with the 
corresponding lower half, the four diagrams (c), (f-h) will emerge. 
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Figure 6. Transition amplitudes to first order as cut Keldysh diagrams 



4-2. Structure of higher- order contributions 

Already the above GW-like approximation gives at least a qualitative description of how the 
leading plasmon satellite is modified by interference. However, although GW theory is usually 
quite accurate for quasi particle positions it gives a somewhat poor representation of the satellite 
shape. In the spectral function A(e) there is only one satellite which is too broad and with a 
mean energy too far away from the quasi particle. By adding the first vertex diagram and self- 
consistency diagrams one improves the satellite in the case of core electrons, and model studies 
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Figure 7. Transition amplitudes contributing to the second plasmon satellite 
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Figure 8. Transition amplitudes modifying the first plasmon satellite 



indicate that also the valence electron spectrum is improved. A corresponding approximation 
for photoemission will involve at least two screened interactions. 

Starting with second order theory, a large number of diagrams occurs which most probably 
gives very small contributions. The screening of the photons can now in principle interfere 
with the excitations shaken up by the photoelectron. If such contributions are neglected, the 
coupling to the photon field enters solely in form of a time-ordered screened optical vertex A as 
in Eq. ([270. Except possibly very close to threshold, the important processes are those where 
photoelectron, once created, remains until it is detected. In terms of diagrams this corresponds 
to adding screened interaction lines to the basic independent-electron triangle in Fig. ^ The 
leading-order amplitudes which contribute to the second plasmon satellite are illustrated in 
Fig. d The dangling interaction lines can be combined in two different ways when complete 
photoemission diagrams are formed. As a result, the amplitudes in Fig. produce 18 complete 
diagrams. The contribution from these 18 diagrams can be written 



Ji 2 ' 2) 



2*£ 



nfiu 



cLdidu>2 



l=a,b,c 



x a 2 n (e p + u)i +lo 2 -u)b 2 Jui)bl(u 2 ), 



(39) 



where 



M 2 X 



(Xp |AG c (e p - Ld)v^(u)i)G c {e p + u\ - u)\v u (uj2)u n (e p + lo x + lq 2 - uj)) 

(Xp v^(uji)\G c (e p + uj 1 )AG c (e p + u\ - u)\v u (uj2)u n (e p + u\ + lo 2 - w)) (40) 

{Xp v^{u:i)\G c (e p + to l )v v {u:2)G c {e p + loi + U2)A\v u (u 2 )u n (e p + u\ + uj 2 - u)) 



In addition, there are amplitudes which modify the shape of the first plasmon satellite 
containing one dangling interaction line (Fig. EJ). The latter amplitudes are to be added to 
the first-order amplitudes which are the squared and summed. The diagrams which result from 
pairing off with the amplitudes Mi i in Fig. El are of second order in W . In this way 24 complete 
Keldysh diagrams are obtained, and they give a contribution of the form 

r°° r i * 

= A^Y, M K^^e^ + ^-^N, (41) 



where Ml 2 is the sum of the six amplitudes M\ 2 / in Fig. |HJ The sum jj 1 ' 1 ) + J^ 2a) and 

thereby the complete second-order current Jp 1 ' 1 ^ + Jp 1 ' 2 ^ + Jp 2 ' 2 ^ is not positive-definite. If 
we add the diagrams where the second-order amplitudes to the first boson satellite are paired 
against themselves a manifestly positive contribution to the first boson satellite is obtained, 

r 00 2 
= 2vr^/ d^M^ + M^ aKep + un-^bfax). (42) 

(1 2) (2 2) 

The sum J p + Jp includes all second-order and part of the third-order contributions to the 
photoelectron current. 

The above procedure can evidently be extended to any order. The basic idea is to sum 
transition amplitudes to a given order, and then form all diagrams containing these amplitudes 
by pairing as above. For a given order n the current is the sum Jp 1 '^ + Jp 2 ' n ^ + • • • + Jp n ' n ^ 
corresponding to 1, 2, . . . , n boson losses. Each partial current Jp' n ^ will involve amplitudes with 

r dangling and up to n — r time-ordered interaction lines. The partial current Jp' n ^ is obtained 
by pairing these amplitudes in all possible ways. The partial current will include all diagrams up 
to order re, and part of the contribution of order between n and 2n — r and it will be manifestly 
positive. If one could continue indefinitely, eventually all diagrams will be included. The method 
of expanding transition amplitudes rather than spectral densities is also applicable for obtaining 
approximations to the one-electron spectral function and to the dynamical structure factor. 



4-3. Model calculations of plasmon and electron-hole losses 

I here describe model calculations in the simplest non-trivial approximation in which we keep 
the four diagrams (c) and (/ — h) in Fig (JIJ. 

In the expressions for matrix elements, G c one can usually approximate it by G r in Eq. 1)36(1 
and by a G a in Eq. (|34|) . In our model calculations we have further approximated G r l a by 
l/(ho ± iT) with a damping from the on-shell value of SjS. The matrix elements can then be 
evaluated by propagating the orbitals in space under influence of a non-Hermitian Hamiltonian 
which greatly simplifies the calculations. The system was approximated by a a semi-infinite 
jellium perturbed by a weak lattice pseudo potential. (Without a lattice potential only surface 
photoemission can occur because of momentum conservation.) The dielectric function and thus 
B(co) was approximated by a semi-analytical result by Bechstedt et al. 42 . 

In Fig. we show results for plasmon losses with an excitation energy 7.5 Ry above 
threshold. These and the following results correspond to aluminium, and the lattice G vector, 
which acts as momentum source, has a parallel component of 0.7 a.u. and a component of 
2.6 a.u. perpendicular to the surface. The extrinsic losses dominate, and the total loss is 
somewhat smaller than the sum of the extrinsic and intrinsic ones. However, we have observed 
that interference depend on the direction of the photocurrent even at 7.5 Ry. The two loss 
mechanisms are therefore not completely decoupled even at 100 eV. In Fig. (|10|) we have 




Figure 9. Photocurrent as function of photoelectron energy in Ry. The zero level corresponds 
to the bottom of the valence band. The spectrum is excited 7.5 Ry above threshold, and the 
photoelectron mean free path is approximated by a constant corresponding to T p = 0.3 Ry. 




Figure 10. Photocurrent excited 4.5 Ry above threshold. T p = 0.3 Ry. Notations as in Fig. 



lowered the excitation energy to 4.5 Ry. In this case the interference is much stronger, and the 
extrinsic and intrinsic parts interfere destructively. 

Finally one may ask how the quasi particle line shape is modified by the transport and 
interference effects. In our calculations we have obtained almost negligible effects. As an 
example, I show the quasi particle region at an excitation energy of 2.5 Ry (Fig. The 
broadening that is seen without any particle-hole effects comes from the finite mean free path 
which smears the momentum selection rules normal to the surface. 

As remarked above, the first-order GW theory gives a rather crude representation of the shape 
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Figure 11. Broadening of the quasi particle peak at normal emission. Ui = 2.5 Ry. For 
notations see Fig. |5] 



of the plasmon satellites. The description will be improved by going to second order in transition 
amplitudes as described in Sec. 14.21 This higher-order approximation is still simple enough to 
allow for numerical evaluations. Within the approximations used above for the first-order results 
the second-order amplitudes can again be obtained by propagating one-electron orbitals which 
makes the calculation manageable. As the excitation energy is further increased, more and more 
plasmon satellites will result in experimental spectra. At sufficiently high energies, the energy 
losses can be described by semi-classical approaches in which the photoelectron is considered 
as a classical particle which moves along a certain trajectory |431 135] or by an integral equation 
which approximates the losses at high energies j!]. We expect that explicit quantum- mechanical 
calculations up to second order in W will cover the energy window up to the limit where these 
simpler methods are applicable. Extensions to second order along the lines given here are on 
the way and will be presented elsewhere. 

5. Concluding remarks 

In this paper I have discussed how techniques from non-equilibrium Green's function theory 
may be used to develop useful approximations as well as for establishing general results such 
as the structure of the no-loss quasi particle part of the spectrum. Although I have discussed 
energy losses in some detail in the second part, I would like to emphasis that it is the no- loss 
quasi particle part that is usually of primary interest in experimental studies. The current 
computational schemes in which the photoelectron moves in a non-Hermitian optical potential 
is essentially correct. The necessary ingredients are here sufficiently accurate approximations for 
the self-energy including its imaginary part. It may be necessary to account for its non-locality 
at least in some approximate way. The proper evaluation of matrix elements the screened optical 
field is another area where improvements may be necessary. Especially in spectral regions where 
the photoelectron self-energy is rapidly varying it may also be necessary to account for the vertex 
corrections discussed in Sec. |3 

In order to treat energy losses below the main quasi-particle peak, a regrouping of diagrams 
in many-body perturbation theory has been proposed. In each order the partial sums have a 
Golden-Rule-like form which guarantees positive spectral intensities. This way of regrouping 



No hole broadening 




may be useful also in other problems involving spectral intensities. One such problem is how to 
improve the one-electron spectral function beyond GW theory. 

Some new model results for valence-electron photoemission have been presented, and more 
realistic results are on the way. The model results confirm the expectations that interference 
effects modify the shape and strength of plasmon satellites also quite far away from the excitation 
threshold, whereas the quasi-particle part is little modified. Based on our ongoing work we intend 
to compare fully quantum-mechanical results to semi-classical descriptions of the photoelectron 
transport away from threshold. 
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